Evolution of sub-horizon cold dark matter perturbations 
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We investigate the evolution of sub-horizon cold dark matter perturbation in the dark energy 
dominated era of the Universe. By generalizing the Meszaros equation to be valid for an arbitrary 
equation of state parameter, we derive the ui-Meszaros equation. Its solutions determine the evolu- 
tion of the cold dark matter perturbation by neglecting dark energy perturbations. Our analytical 
results provide a qualitative understanding of this evolution. 



I. INTRODUCTION 

The cosmological principle is the cornerstone of mod- 
ern cosmology, it assumes that the Universe is homoge- 
neous and isotropic. The cosmological principle is valid 
for studying the Universe as a whole but we know that it 
does not hold perfectly. The nearby Universe is highly in- 
homogeneous, it is made up of stars, planets and galaxies 
rather than a smoothly distributed fluid. 

One of the merits of inflation is that it provides a 
theory of inhomogeneities in the Universe, which may 
explain the observed structures. These inhomogeneities 
arise from the quantum fluctuations in the inflaton field 
about its vacuum state, in other words, by the vacuum 
fluctuation 

The vacuum fluctuations generate a primeval density 
perturbation, that is Gaussian and adiabatic, and whose 
spectral index is close to 1. Such a perturbation was 
regarded, even before the study of inflation, as a viable 
candidate for the origin of the large-scale structure and 
the cosmic microwave background (CMB), which was not 
been observed at that time Q. To understand the evo- 
lution of the primeval perturbation to the present, we 
need to know the nature and amount of the non-baryonic 
dark matter, as well as the value of the cosmological con- 
stant. The simplest possibility is to have zero cosmologi- 
cal constant and cold dark matter giving critical density. 
That result can fairly be said to be the simplest plausi- 
ble model for the origin of large-scale structure and the 
CMB anisotropy and this is called the cold dark matter 
model (CDM) 

Setting up a system of equations to be solved and the 
initial conditions for the perturbations, we can calcu- 
late the inhomogeneities and anisotropies in the Universe. 
The dark matter perturbations are, in principle, coupled 
to all other perturbations. In practice, however, per- 
turbations to the dark matter depend very little on the 
details of the radiation perturbations. Dark matter is 
affected by radiation only indirectly, through the gravi- 
tational potentials. At late times, when the Universe is 
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dominated by matter, these potentials are independent of 
radiation. At early times the gravitational potentials are 
dominated by radiation. The radiation perturbations are 
relatively simple, so that all multipole moments beyond 
the monopole and dipole can be neglected. The converse 
is not true: to treat the anisotropies properly we will 
need to know how the matter perturbations behave. 

One of the earliest quantitative investigations of the 
importance of the gravitational effects of CDM on the 
baryonic fluctuations which give rise to galaxies was pre- 
sented by Meszaros in These papers discussed the 
evolution of the density contrast in what came to be 
known as cold dark matter fluctuations in the cosmo- 
logical radiation background as it enters the horizon and 
goes through the matter-radiation equilibrium and the 
recombination epochs, and its effect on the formation of 
galaxies and large scale structure mediated by such cold 
dark matter. This is referred to in the literature as the 
"Meszaros effect" , which has become incorporated in all 
subsequent CDM calculations of structure evolution, de- 
termining the so-called transfer function of the original 
fluctuation spectrum. Large scale structure observations, 
modeling and simulations are now a major activity in cos- 
mology, but the actual nature of dark matter still remains 
elusive. 

In this work, we want to study the evolution of the 
dark matter perturbations in the dark matter-dark en- 
ergy era. In Section [II] we review the Meszaros equation 
and its well known solutions. In Section IlIII we obtain a 
Meszaros-like equation: the w-Meszaros equation, by us- 
ing a similar formalism used to obtain the original equa- 
tion. In Section IIVI we find all solutions of this equation 
and in the final Section [V] we discuss and analyze our 
results. 



II. MESZAROS EQUATION 

The radiation pressure causes the gravitational poten- 
tials to decay as modes enter the horizon during the ra- 
diation era. The pressure suppresses any growth in the 
radiation perturbations. In contrast, the matter pertur- 
bations grow logarithmically. In the radiation era, the 
potential is determined by radiation but eventually the 
growth in the matter perturbations offsets that there is 
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more radiation than matter. The dark matter perturba- 
tions (5 C ) become larger than the radiation perturbations 
(S r ) even if p c is smaller than p r . Once this happens, 
the gravitational potential and the dark matter pertur- 
bations evolve together and decouple from radiation. 

The perturbation equations for dark matter in sub- 
horizon scales are: 
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where the derivative is with respect to the proper time. 
Note that in the Poisson Eq. ([3]) we are not taking into 
account the radiation perturbations, even when 4> is de- 
termined by them, by this time S c > S r and as we pointed 
before, <fr and S c evolve together. 

We want to follow the perturbations through the epoch 
of equality, so is better to define a new temporal vari- 
able 0]: 
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and in terms of this new variable the equations become, 
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where the prime denotes differentiation with respect to 
y. By turning this first order system into a second-order 
differential equation for S c we obtain the Meszaros equa- 
tion which governs the evolution of sub-horizon cold dark 
matter perturbations once radiation perturbations have 
become negligible: 
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This equation has two solutions given by Meszaros 
The growing solution is, 
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and the decaying solution is given by 
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Before equality (y <C 1), the growth mode is practically 
frozen, whereas at late times (y ^> 1) the growing solu- 
tion Sc^ scales as y and the solution matches the law in 
a matter dominated universe 
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The decaying mode S c falls off as y 
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III. w-MESZAROS EQUATION 

Now we want to study the evolution of the dark matter 
perturbations in the dark matter-dark energy era. We 
are now going to neglect the dark energy perturbations, 
analogously to neglecting the radiation perturbations in 
the Meszaros effect previously. 

Along the lines of ((H) we introduce a new dimensionless 
variable as an evolution parameter instead of t or a(t) 
which we define by 
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where w is the constant equation of state parameter of 
the dominant component of the Universe, p m corresponds 
to cold dark matter p c plus baryons pb and a co is the 
coincidence time: ^^p) ~ 1- It should be noted that now 

y is a decreasing functions with respect to the expansion. 

In terms of this variable equations dU, ([2]) and ([3]) 
become 
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We are going to differentiate Eq. ([T 
into a second-order equation 
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Using the velocity equation Eq. (1141) we can eliminate 
v' m and the continuity equation (|13[) to eliminate the v m - 
term and from 
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dy \aHy 
the equation becomes 
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By using the Poisson equation (|15p we can finally elimi- 
nate the gravitational potential 4> and arrive at 
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We refer to this as the w-Meszaros equation for dark 
matter-dark energy, which governs the evolution of sub- 
horizon cold dark matter perturbations if we neglect dark 
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energy perturbations. This equation tells us that in the 
dark matter-dark energy era, the dark matter perturba- 
tions will decrease as the dark energy starts to dominate. 
If we set w = 1/3, for a radiation dominated Universe, 
we recover the well-known Meszaros equation. 



IV. SOLVING THE w-MESZAROS EQUATION 

Let us introduce the new independent variable z — —y, 
then (|19p can be written in the form 
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which is of the form of the hypergeometric differential 
equation 

z(l - z)44 + [c- (a + b+l)z]^f- - abS = 0. (21) 
dz z dz 

This differential equation has two linearly independent 
solutions which both involve the hypergeometric function 

W\ h 1 T(c) ^ T(a + n)T(b + n) z n 

F[a,b,c;z\ = > — — — r -, (22 
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where T is the Gamma function. The exact form of the 
solution depends on the numerical values of the constants 
a, b, c, see [§, Chapter 15]. 

For the w-Meszaros equation the three constants a, b 
and c can be expressed in terms of w only, and are given 
by 
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c= — (27w 3 - 9w 2 + 2) 
6w 

d = 2916w 6 - 972w 5 - 567w 4 + 324w 3 - 24w + 28. 

Now, if we choose w = 1/3 then a = 1, b = —3/2, c = 0, 
in which case the general solution involving the hyperge- 
ometric function reduces to Eqs. ([9]) and (|10[) . In princi- 
ple one can now study solutions of (|19p numerically for 
arbitrary w. 

However, for some special values of the equation of 
state parameter, we can find explicit solutions. We found 
that for the two choices w = —1/3 and w = 1/6 this 
is possible, and we could not find other values in the 
parameter range — 1 < w < 1, including w = — 1. For 
w = —1/3 the solution is 
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which only contains elementary functions. The c\ and c 2 
part in this solution are the growing and decaying modes 
with respect to y, respectively. This solution is similar to 
the well-known w — 1/3 solution, see Eq. © and ([TP")) . 
The solution for w — 1 /6 is given by 

S(y) = ci(3 + 144y + 256y 2 )(l + l/y) 7 ^ 8 

+ c 2 (l + l/y) 7 / 16 Ql /S (l + 2y), (24) 

where is the generalized Legendre polynomial and we 
are interested in the real part of this solution. As for the 
previous solution, the c± and ci part are the growing and 
decaying modes, respectively. 

When the equation of state parameter is w = —2/3, we 
have a physically particularly interesting solution since 
this models one form of dark energy. One cannot write 
the solution in terms of elementary functions. In terms 
of the hypergeometric function, we can write 
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where the constants c\ and ci are specified by initial or 
boundary conditions. Note that the independent variable 
y and the scale factor a are related by y = 1/a 2 . 



V. ANALYSIS AND DISCUSSION 

Despite the fact that solutions of the iu-Meszaros equa- 
tion (fl9|) are rather complicated, the three constants a, 
b and c depend only on the value of the equation of state 
parameter w. One can easily analyze these solutions nu- 
merically and for instance plot the late time evolution of 
the dark matter perturbations. The left panel of Fig. Q] 
shows the evolution of the perturbations 5(a) of Eq. (|25|) 
from the dark matter dominated era to the distant fu- 
ture, the dark energy dominated era, for w = —2/3. The 
a = 1 line corresponds to today. As expected, the dark 
matter perturbations start to decrease after the time of 
the coincidence. 

This solution is characterized by two decaying modes 
which behave rather differently. One part (the c 2 part) 
of solution (|25p decays for small values of a and then 
approaches a constant at late times. The other part (the 
ci part) of the solution stays constant at early times and 
rapidly decrease at some point in the near future. This 
part dominates the late time behavior of the solution 
which is clearly illustrated in Fig. [T] (left panel) . When 
we further decrease the value of w, we find the same 
qualitative picture, the right panel in Fig. [T] shows the 
solutions for w = — 1. 

As we pointed before the dark matter perturbations 
start to decrease after the time of the coincidence and will 
no longer evolve as S c ~ a. At some time in the future, 
when the Universe gets dominated by dark energy they 
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FIG. 1: The left panel shows the numerical evolution for w — —2/3, the right panel shows the same scenario for w = — 1. We 
have chosen the boundary conditions S(y — 10~ 4 ) = 0.2 and 5(y = 1) = 1. The solid line represents the complete solution, 
the dotted and dashed lines correspond to the ci and C2 parts of the solutions respectively. In both plots the dotted lines have 
been moved to 0.1 by hand to make them visible. 



will become negligible and the dark energy perturbations 
will dominate. This time will depend on the equation of 
state of the dark energy, i.e. will depend on the model of 
dark energy or modified gravity considered to solve the 
equation. Therefore, the w-Meszaros equation can be 
used to distinguish between different cosmological mod- 
els. For practical purposes, and to get some predictions 
out, one would have to use current data and match the 
solutions in order to find a transfer function by following 
well-known methods applied to obtain the growth func- 
tion @. 

It should also be emphasized that one can, in principle, 
use the w-Meszaros equation to study phantom dark en- 
ergy models, where w < — 1, see for example [Iol - [l3j . The 
treatment of perturbations in such models can be quite 
involved, however, our equations are well defined for all 
w. We suggest to use this equation to gain a preliminary 
understanding of phantom dark energy perturbations in 
the late time universe. 

Another interesting idea we would like to point out is 
the derivation of a Meszaros like equation in situations 
where couplings between dark matter and dark energy 



are allowed. This would, as before, allow for studying a 
large class of models qualitatively, see also 

Finally, solving the w-Meszaros equation can give us 
an idea of when the dark energy perturbations become 
important and at what time (size of the universe) they 
have to be taken into account in calculations for an spe- 
cific model. Note that in this work we are not considering 
models with early time dark energy that can lead to some 
early time presence of dark energy perturbations. There- 
fore, it would be very interesting to study the evolution 
of dark energy perturbations using the method described 
here, however, this is beyond the scope of this brief re- 
port. 
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